Spin liquid correlations, anisotropic exchange and symmetry breaking in Tb 2 Ti 2 07 
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We have studied the low energy spin dynamics between 4.6 K and 0.07 K in a Tb2Ti207 single 
crystal sample, by means of inelastic neutron scattering experiments. The spectra consist in a dual 
response, with a static and an inelastic contribution, showing striking Q-dependences. We propose 
an interpretation involving an anisotropic exchange interaction in combination with a breaking of 
the threefold symmetry at the rare earth site. Simulations of the Q-dependent scattering in the 
Random Phase Approximation account well for the inelastic response. 



Spin liquids now attract considerable attention in mod- 
ern condensed matter physics [U • In a classical pic- 
ture, the localized magnetic moments in such coopera- 
tive paramagnets keep fluctuating in a correlated man- 
ner, failing to develop long range order down to very 
low temperature. From a quantum point of view, a spin 
liquid ground state can be described by entangled spin 
wavefunctions and supports exotic fractionalized excita- 
tions also called spinons. Geometrically frustrated mag- 
nets are good candidates in the pursuit of such disordered 
quantum ground states [3] and one of the celebrated ex- 
amples is the Tb2Ti207 pyrochlore. It remains in a spin 
liquid state, with short range correlated fluctuating mo- 
ments, down to a temperature as low as 20 mK [4[ |5]. 
Since 1999, it has been the subject of many theoretical 
as well as experimental works, and the origin of its spin 
liquid ground state is still puzzling. 

Tb2Ti 2 07 belongs to the same family as the Ho2Ti 2 07 
and Dy2Ti207 spin ices, characterized by an Ising 
anisotropy along local < 111 > axes [6 j. However, the 
Tb 3+ crystal electric field (CEF) with trigonal symmetry 
[7J [8] has a much lower energy gap between the ground 
state doublet and the first excited doublet than in spin 
ices (18 K instead of 200 to 300 K). It was suggested that, 
unlike in spin ices, this gap is small enough to allow 
admixture of excited crystal field states, that produces 
an effective ferromagnetic contribution which competes 
with the original antiferromagnetic interactions, and that 
moves Tb2Ti2 07 towards a "quantum spin-ice" regime 
[9] [10] . More recent general descriptions of pyrochlores 
introduce a minimal Hamiltonian, based on symmetry 
grounds, for pseudospins 1/2 (the subspace spanned by 
the ground doublet states \ip±)) [TTHT3] . These models 
involve an Ising exchange constant J zz responsible for 
the spin-ice behavior, as well as three "quantum" terms 
J±i Jz± an d J±± that lift the macroscopic degeneracy of 
the spin-ice manifold and stabilize a so called Coulomb 
phase or U(l) spin liquid phase, describable by an emer- 
gent U(l) gauge field. For large quantum couplings, con- 
ventional phases are stabilized against the spin liquid. 
Interestingly, in the particular case of non-Kramers ions 
(like Tb 3+ ), these ph ases are characterized by ordering 



of the 4/ quadrupoles [TJJ H5j . breaking spontaneously 
the threefold symmetry of the crystal field. 

Recently, we proposed a somehow more phenomeno- 
logical route to point out the relevance of such a sym- 
metry breaking. Indeed, inelastic neutron scattering ex- 
periments have evidenced low energy spin fluctuations 
[8j [16j [17], which, because of general properties of non- 
Kramers ions, cannot be explained in the nominal trig- 
onal CEF, but can be quite naturally accounted for as- 
suming a breaking of the threefold symmetry [T8H20] . 
The same conclusion holds in the case of the ordered 
spin ice parent compound Tb2Sn2 07 [21 , where simi- 
lar, although better defined, strong low energy fluctua- 
tions are also reported [TTJ [22j [23] . Experimentally, such 
a symmetry lowering could be due to a tetragonal dis- 
tortion precursor to a T ~ Jahn- Teller transition. In 
Tb2Ti207, this is supported by thermodynamic measure- 
ments [24, 25 , Raman scattering [26], by some studies of 
the thermal evolution of the elastic constants [27-29 , as 
well as X-ray measurements [30] , but the existence of this 
distortion is still debated in literature [3U [32] . Recently, 
motivated by our work [19] . there have been several at- 
tempts to determine the characteristics of the low energy 
spin dynamics and especially decide if the fluctuations 
are quasi-elastic or inelastic [32j [33], but no consensus 
was obtained. 

In this work, we report new high accuracy and high 
resolution INS experiments performed on a Tb2Ti2C>7 
single crystal down to 0.07 K. We unambiguously ob- 
serve the controversial inelastic excitation. More pre- 
cisely, we found that the low energy response is dual, i.e. 
it is the sum of a rather strong elastic signal and of an 
overdamped inelastic one with striking Q-dependences. 
By modeling the dynamic susceptibility in the Random 
Phase Approximation (RPA) and using the anisotropic 
exchange tensor determined previously [19 , we find that 
the hypothesis of a breaking of the threefold symmetry 
reproduces the Q-dependence of the inelastic response 
satisfactorily. 

INS experiments have been performed on the triple 
axis spectrometer 4F2, installed on the cold neutron 
source at the LLB-Orphee (Saclay, France) neutron facil- 
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FIG. 1: (color on line) a: Neutron intensity (normalized 
to monitor) as a function of energy transfer at T=0.07K for 
Q=(002). The data were recorded with a final wavevector 
k/=1.2A _1 , yielding an energy resolution (FWHM) Ao = 
0.07 meV. The lines show the elastic and inelastic contribu- 
tions; the red one shows the total fitted intensity in the case of 
a damped harmonic oscillator (DHO) for the inelastic contri- 
bution (see text). Insert: thermal evolution of the low energy 
scattering, b and c: Elastic magnetic scattering maps taken 
respectively at 0.07 K and 4K. A Q-independent background 
has been subtracted from the raw data in order to remove the 
incoherent scattering. It was estimated from the temperature 
dependence of the total elastic scattering. The color scale is 
identical for the two temperatures. The red box corresponds 
to the points actually measured, the remaining has been de- 
duced by symmetry. 



ity (see the supplemental material for detailed informa- 
tion). Four Tb2Ti207 single crystals, of total mass 11 g, 
synthesized by the floating zone method, were co-aligned 
in the (hhl) scattering plane. We collected a series of 
energy scans at various wavevectors along the high sym- 
metry directions (/i/i0), (00/) and (hhh) for different tem- 
peratures ranging from 0.07 K up to 4 K. We also mapped 
out the intensity at different constant energy transfers as 
a function of wavevector in the (hhl) plane. Our data 
provide compelling evidence for a strong low energy re- 
sponse, well below the first CEF transition (at about 
1.5 meV). The high energy resolution Ao=0.07meV (full 
width at half maximum FWHM) allowed us to separate 



an elastic contribution Iq(Q) from a broad inelastic con- 
tribution ii(Q,cj), with a maximum around 0.15 meV. 
This is illustrated in FigjTJ showing the neutron intensity 
at Q=(002) in the low energy range —0.1 < uo < 0.4 meV. 
The two contributions are well separated at 0.07 K, where 
a clear shoulder is observed at finite energy transfer. 
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FIG. 2: (color on line) Inelastic scattering maps at 0.07 K for 
an energy transfer: a: 0.08 meV , b: 0.11 meV, c: 0.15 meV 
and d: 0.20 meV. These energies are indicated by arrows in 
figure [TJl. 

This dual response was modeled by I(Q,lj) = Io(Q) + 
I\(Q,uj). Note that because of the experimental res- 
olution, Iq(Q) either corresponds to pure static or to 
slow fluctuations with characteristic frequency lower than 
v ^ 4 GHz. Like in frustrated ferromagnets, dif- 

ferent forms have been used for I\(Q,uj): the double 
Lorentzian (DL) profiles and the damped harmonic oscil- 
lator (DHO), with characteristic energy oj\ and damping 
(HWHM) Ti [34 . After convolution with the experi- 
mental resolution function, we found that the DHO pro- 
file yields the best agreement with the data. Fitting the 
energy scans recorded along the high symmetry direction 
shows almost no change of uo\ in Q-space: uo\ ~ 0.20 meV, 
while the damping Y\ is large and has the same magni- 
tude: Ti « uj\. Maps of Iq(Q) measured at 0.07 K and 
4 K, shown in Figjl] reveal a ground state with strongly 
anisotropic short range correlations. At 0.07 K, we ob- 
serve a complex pattern resembling an array of " butter- 
flies" pinned at (002) and, with a smaller intensity, at 
(220)-type positions. Superimposed on that structure, 
the map shows lobes centered around (1/2, 1/2, 1/2)- 
type positions, elongated along the (hhh) directions, and 
separated by a pinch-point at (111). A series of inelastic 
Q-space maps, taken at 0.07 K at different energy trans- 
fers, are shown in Figj2] As the energy transfer increases, 
the pinch-point at Q=(lll) is progressively "filled" while 
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the intensity at (l/2,l/2,l/2)-type positions decreases. 
With increasing temperature up to 4K, the intensity of 
the two signals progressively weakens and the inelastic 
signal progressively merges into the elastic one. The lobes 
centered at (1/2,1/2,1/2) persist up to 0.5 K, but become 
hardly visible at 4K, with only weak maxima at (002) 
and (220)-type positions. 

We proceed now to the analysis of the inelastic maps in 
the (hhl) plane of the reciprocal space. As mentioned in 
the introduction, the hypothesis of a breaking of the trig- 
onal symmetry, possibly due to a precursor Jahn- Teller 
distortion developing at low temperature, has been put 
forward to explain the peculiar properties of Tb2Ti 2 07 
[H1[19]. The relevant CEF Hamiltonian writes: 

HCEF = Htrig + ^[2JJ + Jl + V2(J X J Z + J Z J X % (1) 

where Htrig represents the nominal trigonal crystal field 
and the second term a small tetragonal distortion along 
a cubic [001] axis (J is the total angular momentum; 
for Tb 3+ , J=6 and the Lande factor is gj = 3/2). The 
breaking of the threefold symmetry results in a degener- 
acy lifting and in a mixing of the two states of the ground 
doublet l^i), and stabilizes a CEF singlet state described 
by the symmetric wavefunction |^ s ) = ~^[| < + ) + |-0_)] 
P~8j [19] . The first excited state is then the antisymmet- 
ric combination: \i/; a ) = ^[|^ + ) — 1^-)], separated from 
\ip s ) by a quantity S proportional to the amplitude of 
the distortion. In the neutron spectra, this gives rise 
to an excitation at an energy close to <S, unveiling the 
transition from \ijj s ) to \^ a )- Its intensity is propor- 
tional to ^2 a=x y z \(ip s \J a \^a}\ 2 • A straightforward cal- 
culation shows that, whereas the general properties of 
non-Kramers ions impose (ip+\J\ip_) = 0, the entangle- 
ment of \ip±) yields a sizeable (^a\J\^s) matrix element 
and thus a large cross section. This accounts for the 
presence of the low energy inelastic line in Figjl] a. In 
the absence of entanglement, its intensity would be van- 
ishingly small [23] . 

We then introduce the exchange/dipolar Hamiltonian 
widely accepted for pyrochlores [7], with an anisotropic 
exchange tensor J ex taken to be diagonal in the (u, v, w) 
frame linked with a Tb-Tb bond along w [35 . For in- 
stance, for the bond along [110], this frame is defined 
by: £=(0,0,1), v= ^(1,-1,0) and w = ^(1,1,0). The 
link between this exchange tensor and that derived in 
Ref.[T3] is discussed in the Supplemental material. For 
not too large antiferromagnetic J ex , depending on the 
value of the J ex /5 ratio, the mean field ground state of 
the model is either an ordered spin ice phase or a sin- 
glet yielding no magnetic order, which is quite similar to 
Bleaney's result for non-Kramers rare earth ions with a 
singlet CEF ground state [36] . 

The set of parameters {Dq and the exchange ten- 
sor components) that stabilizes \ip s ) as the ground state 
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FIG. 3: (color on line) Inelastic scattering maps in Tb2Ti 2 07. 
Left column: T=4.6K and energy transfer 0.15 meV: a: Ex- 
perimental data; b and c: Simulations respectively with a 
tetragonal distortion Dq — 0.25 K and without. The simula- 
tion in c has been multiplied by 7000 to obtain the same inten- 
sity as in b. Right column: T=0.07K and energy transfer 
0.08 meV: d: Experimental data; e and f: Simulations in the 
presence of a tetragonal distortion Dq = 0.25 K, with respec- 
tively an anisotropic and an isotropic exchange tensor. From 
these latter maps, it is clear that the Q-dependence of the 
scattering strongly depends on the anisotropy, and that it is 
not reproduced with isotropic exchange. 



has been determined previously [T9] : Dq=0.25K, J u = 
-0.07 K, J v = -0.2 K and J w = -0.1K (a negative 
sign for Ji means an antiferromagnetic coupling). With 
this Dq value, a low energy inelastic line is expected at 
8 ~ (ji ~ 0.22 meV in the inelastic neutron data. We 
then calculated the corresponding cross section in the 
RPA approximation [37-39J as described in the Supple- 
mental material. A finite width (HWHM) Fi=0.20meV 
was introduced in the single site susceptibility for the 
lower transition \ip 3 ) <o> |^a), according to the DHO fit. 
The widths for the other transitions were set at a smaller 
value, with no influence on the result. 

Figure[3]compares the experimental inelastic scattering 
[3^i at 4.6 K for an energy transfer 0.15 me V with simula- 
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tions obtained with ^p) and without ^p) the tetragonal 
distortion. The simulated map with distortion is much 
closer to the experimental data, reproducing the two in- 
tense spots at (002) and the rhomb-shaped lesser inten- 
sity scattering centered at (220). The map simulated 
without distortion, i.e. assuming a degenerate doublet as 
ground state, is quite different, but it bears resemblance 
with the diffuse scattering pattern observed at 9 K [5j [37] . 
This points to the fact that the symmetry breaking is al- 
ready present at 4 K, in agreement with the Raman data 
[26] . At the base temperature (70 mK), the simulated 
inelastic map [3^, for an energy transfer 0.08 meV, also re- 
produces correctly the structure in Q-space reproduced 
for convenience in [3]d: the observed butterfly- shaped 
structures with centers at (002), the high intensity lines 
along (111) and the pinch-points at (lll)-type positions 
are present. Finally, we examine the influence of ex- 
change anisotropy. Figure [3]f represents a simulation at 
0.07 K performed in the same conditions as[3^, but for an 
isotropic exchange constant J = — 0.04K [19]. Its struc- 
ture does not show the "butterfly-shaped" pattern ob- 
served in the data. Our calculations actually reveal that 
such a pattern arises only for an anisotropic exchange 
tensor close to that derived for Tb 2 Ti 2 7 (J u = -0.07 K, 
X = -0.2K, J W = -0.1K). 

These RPA simulations thus point to the relevance of 
the symmetry breaking, leading to a sizeable inelastic 
scattering, while the anisotropy of the exchange coupling 
does reproduce quite well the Q-dependence of the in- 
elastic response. We thus have shown that at this mean 
field level the model captures important features relevant 
to the physics of Tb2Ti207. It is however inadequate to 
handle other complex features. Especially, the elastic 
component, given by the Curie term in the expression of 
the single site susceptibility (see Supplemental Material), 
is zero since the CEF singlets are non magnetic. This is a 
limitation of the mean field approach, but we notice that 
experimentally, the Q dependence of the elastic signal is 
similar to that of the inelastic one, suggesting that the 
same set of interactions is at play. 

To go beyond the present model, different approaches 
could be suggested. At the level of a single tetrahedron, 
exact diagonalization of the Tb spin system, currently 
under way, could possibly capture part of the elastic scat- 
tering. Coupling the moments by defects, such as a dis- 
tribution of distortions for instance, could yield static 
correlations and the exotic transition seen in suscepti- 
bility and specific heat[25j 00] [41]. Other approaches 
suggest that tunneling processes between different spin- 
ice configurations also give rise to quasi-elastic scatter- 
ing [IT] . Finally, we note that recent theoretical models 
for non-Kramers ions predict new phases characterized 
by a spontaneous breaking of the threefold symmetry of 
the crystal field [14] [15] as well as an ordering of the 4/ 
quadrupoles. Since a coupling between lattice degrees of 
freedom and such quadrupoles is quite natural, this could 
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FIG. 4: Simulated diffuse SF (left) and NSF (right) scattering 
maps at 0.07 K in Tb 2 Ti 2 7 . 



provide a new basis for our interpretation. 

After this paper was submitted, we became aware of 
an independent neutron scattering study by Fennell et 
al [42] . It was performed in different experimental con- 
ditions, being an energy integrated experiment, whereas 
our measurements are performed with energy analysis of 
the outcoming neutrons, thus yielding a direct access to 
the spin fluctuation spectrum. Fennell et al also report 
on polarized neutron scattering experiments, which allow 
to distinguish the correlations of spin components along 
the vertical axis [1,-1,0], in the Non-Spin-Flip (NSF) 
channel, and correlations of spin components perpendic- 
ular to Q within the scattering plane in the Spin-Flip 
(SF) channel. Figure [4] shows, for both channels, energy- 
integrated maps at 0.07 K calculated in the framework of 
our model, that capture the experimental data (Fig. 2 
A and F of ref.[42]) quite well. The hypothesis of a sym- 
metry breaking and anisotropic exchange are thus com- 
pulsory to understand the inelastic response, as we have 
shown above, but also the diffuse scattering observed in 
ref. [42]. 

In conclusion, our INS experiments in Tb2Ti207 per- 
formed down to 0.07 K show unconventional spin dynam- 
ics with peculiar features in Q-space, such as pinch-points 
and "butterfly-shaped" patterns in the (hhl) plane, very 
different from those observed in classical spin ices. We 
propose an interpretation in terms of a breaking of the 
threefold symmetry of the crystal field at the Tb 3+ sites. 
We present calculations of the inelastic scattering and 
diffuse cross-section of polarized neutrons in the frame of 
the RPA approximation which support this picture and 
confirm the anisotropy of the exchange tensor derived 
previously for Tb2Ti 2 07. 

We would like to acknowledge fruitful discussions with 
P. Holdsworth, M. Gingras, B. Canals, E. Lhotel as well 
as with B. Hennion. 



5 



[1] C. Lhuillier and G. Misguich, in Introduction to frustrated 
magnetism, Eds. C. Lacroix, Ph. Mendels and F. Mila 
(Springer 2011) 

[2] L. Balents, Nature 464, 199 (2010) 

[3] J. S. Gardner, M. J. P. Gingras and J. E. Greedan, Rev. 
Mod. Phys. 82, 53 (2010) 

[4] J. S. Gardner, S. R. Dunsiger, B. D. Gaulin, M. J. P. 
Gingras, J. E. Greedan, R. F. Kiefl, M. D. Lumsden, W. 
A. MacFarlane, N. P. Raju, J. E. Sonier, I. Swainson and 
Z. Tun, Phys. Rev. Lett. 82, 1012 (1999). 

[5] J. S. Gardner, B. D. Gaulin, A. J. Berlinsky, P. Waldron, 
S. R. Dunsiger, N. P. Raju and J. E. Greedan, Phys. Rev. 
B 64, 224416 (2001). 

[6] H. Cao, A. Gukasov, I. Mirebeau, P. Bonville, C. Decorse 
and G. Dhalenne,Phys. Rev. Lett. 103, 056402 (2009) 

[7] M.J. P. Gingras, B.C. den Hertog, M. Faucher, J.S. Gard- 
ner, S.R. Dunsiger, L.J. Chang, B.D. Gaulin, N.P. Raju, 
J.E. Greedan, Phys. Rev. B 62, 6496 (2000) 

[8] I. Mirebeau, P. Bonville, M. Hennion, Phys. Rev. B 76, 
184436 (2007) 

[9] H. R. Molavian, M. J. P. Gingras, and B. Canals, Phys. 
Rev. Lett. 98, 157204 (2007) 
[10] B. C. den Hertog, M. J. P. Gingras, Phys. Rev. Lett. 
84,3430 (2000) 

[11] O. Benton, O. Sikora and N. Shannon, [arXiv: 1204. 1325 
(2012) 

[12] M. Hermele, M. P. A. Fisher and L. Balents, Phys. Rev. 

B 69, 064404 (2004) 
[13] L. Savary and L. Balents, Phys. Rev. Lett. 108, 037202 

(2012) 

[14] Sung Bin Lee et al, arXiv. 1204.2262 

[15] S. Onoda and Y. Tanaka, Phys. Rev. Lett. 105, 047201, 

(2010); S. Onoda and Y. Tanaka, Phys. Rev. B 83, 

094411 (2011). 

[16] Y. Yasui, M. Kanada, M. Ito, H. Harashina, M. Sato, H. 

Okumura, K. Kakurai, and H. Kadowaki, J. Phys. Soc. 

Jpn. 71, 599 (2002) 
[17] K. C. Rule, G. Ehlers, J. S. Gardner, Y.Qiu, E. Moskvin, 

K. Kiefer and S. Gerischer, J. Phys.: Condens. Matter 

21, 486005 (2009). 
[18] P. Bonville, I. Mirebeau, A. Gukasov, S. Petit and J. 

Robert, J. Phys.: Conf. Series 32, 012006 (2011). 
[19] P. Bonville, I. Mirebeau, A. Gukasov, S. Petit, J. Robert, 

Phys. Rev. B 84, 184409 (2011) 
[20] K. C. Rule, P. Bonville, J. Phys. Conference Series 145 

012027 (2009) 

[21] I. Mirebeau, A. Apetrei, J. Rodriguez-Carvajal, P. 
Bonville, A. Forget, D. Colson, V. Glazkov, J. P. Sanchez, 
O. Isnard and E. Suard, Phys. Rev. Lett. 94, 246402 
(2005). 



[22] K. C. Rule, G. Ehlers, J. R. Stewart, A. L. Cornelius, P. 
P. Deen, Y. Qiu, C. R. Wiebe, J. A. Janik, H. D. Zhou, 
D. Antonio, B. W. Woytko, J. P. Ruff, H. A. Dabkowska, 
B. D. Gaulin and J. S. Gardner, Phys. Rev. B 76, 212405 
(2007). 

[23] S. Petit, P. Bonville, I. Mirebeau, H. Mutka, and J. 

Robert, Phys. Rev. B 85, 054428 (2012) 
[24] Y. Chapuis, A.Yaouanc, P. Dalmas de Reotier, C. Marin, 

S.Vanishri, S. H. Curnoe, C. Vaju, and A. Forget, Phys. 

Rev. B 82, 100402(R) (2010) 
[25] A. Yaouanc, P. Dalmas de Reotier, Y. Chapuis, C. Marin, 

S. Vanishri, D. Aoki, B. Fak, L.-P. Regnault, C. Buisson, 

A. Amato, C. Baines, and A. D. Hillier, Phys. Rev. B 84, 

184403 (2011) 

[26] T. T. A. Lummen, LP. Handayani, M.C. Donker,D. 

Fausti, G. Dhalenne, P. Berthet, A. Revcolevschi, and 

P.H.M van Loosdrecht, Phys. Rev. B 77, 214310 (2008) 
[27] L. G. Mamsurova, K. S. Pigal'skii, K. K. Pukhov, JETP 

Lett. 43, 755 (1986) 
[28] Y. Nakanishi et al, Phys. Rev. B 83, 184434 (2011) 
[29] V. V. Klekovkina, A. R. Zakirov, B. Z. Malkin and L. A. 

Kasatkina, J. Phys.: Conf. Series 324, 012036 (2011) 
[30] J. P. C. Ruff, B. D. Gaulin, J. P. Castellan, K. C. Rule, 

J. P. Clancy, J. Rodriguez and H. A. Dabkowska, Phys. 

Rev. Lett. 99 237202 (2007). 
[31] K. Goto, H. Takatsu, T. Taniguchi and H. Kadowaki, J. 

Phys. Soc. Jpn. 81 015001,(2012) 
[32] B. D. Gaulin, J. S. Gardner, P. A. McClarty and M. J. 

P. Gingras, Phys. Rev. B. 84, 140402 (2011). 
[33] H. Takatsu, H. Kadowaki, Taku J. Sato, J. W. Lynn, 

Y. Tabata, T. Yamazaki, and K. Matsuhira, J. Phys.: 

Condens. Matter 24, 052201 (2012) 
[34] B. Hennion, M. Hennion, F. Hippert, A. P. Murani, Phys. 

Rev. B 28, 5365(R) (1983) 
[35] B.Z. Malkin, T. T. A. Lummen, P. H. M. van Loosdrecht, 

G. Dhalenne and A. R. Zakirov, J. Phys.: Condens. Mat- 
ter 22 (2010) 276003 
[36] B. Bleaney, Proc. Roy. Soc. (London) 276a, 19 (1963) 
[37] Y. J. Kao, M. Enjalran, A. Del Maestro, H. R. Molavian 

and M. J. P. Gingras, Phys. Rev. B 68, 172407 (2003). 
[38] J. Jensen and A. R. Mackintosh, Rare earth magnetism 

(Clarendon Press, Oxford, 1991) 
[39] A. G. Del Maestro and M. J. P. Gingras, J. Phys.: Con- 
dens. Matter 16, 3339 (2004) 
[40] E. Lhotel, C. Paulsen, P.D. deReotier, A. Yaouanc, C. 

Marin and S. Vanishri, Phys. Rev. B 86 020410(R) (2012) 
[41] S. Legl, C.Krey, S. R. Dunsiger, H.A. Dabkowska, J. A. 

Rodriguez, G. M. Luke, and C. Pfleiderer, Phys. Rev. 

Lett. 109, 047201 (2012). 
[42] T. Fennell, M. Kenzelmann, B. Roessli, M. K. Haas, and 

R. J. Cava, Phys. Rev. Lett. 109, 017201 (2012). 



